Superconductivity from orbital nematic fluctuations 
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Recent experiments suggest that besides of antiferromagnetic fluctuations nematic fluctuations 
may contribute to the occurrence of superconductivity in iron pnictides. Motivated by this obser- 
vation we study superconductivity from nematic fluctuations in a minimal two-band model. The 
employed band parameters are appropriate for iron pnictides and lead to four pockets for the Fermi 
line. It is shown that low-energy, long-wavelength nematic fluctuations within the pockets give rise 
to strong-coupling superconductivity whereas the large momenta density fluctuations between pock- 
ets are rather irrelevant. The obtained transition temperatures are similar to those typically found 
in the pnictides and are rather robust against repulsive Coulomb interactions. The superconducting 
and nematic states coexist in a large region of the phase diagram. 
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Electronic analogues of nematic liquid crystals attract 
much interest and are observed in a number of con- 
densed matter systems, such as quantum spin systems 
[l[ , two-dimensional electron gases 0, Q , cuprates 0, [B[ , 
bilayer ruthenates [6|, and iron pnictides [7|. In a ne- 
matic state the orientational symmetry is broken whereas 
the other symmetries are retained. Depending on elec- 
tronic degrees of freedom responsible for nematic order, 
we may distinguish between three kinds of nematicity: 
The charge nematicity which is obtained either by par- 
tial melting of charge stripes Q or by a Pomeranchuk 
instability jl,[loj]; the spin nematicity which is driven by 
frustration of spin-spin interactions in quantum spin sys- 
tems or by two competin g sp in-density-wave (SD W) 
order parameters in metals [12| : the orbital nematicity 
due to orbital order caused, for example, by different oc- 
cupations of d yz and d zx orbitals in d-electron systems 

In iron pnictides, nematicity is observed just above 
the SDW phase and is accompanied by a tetragonal- 
orthorhombic structural phase transition [7] . Because of 
the proximity to the SDW phase, one plausible scenario 
is the spin nematicity [15l |16|. On the other hand, an 
occupation difference of d yz and d zx orbitals was directly 
observed in angle-resolved photoemission spectroscopy 
@; E3 • Orbital nematicity is thus equally well a plausi- 
ble scenario 18-22|. Nematic charge fluctuations can be 
observed directly in Raman scattering experiments [23| . 
and were found close to the SDW phase [24j. 

Quite recently, magnetic torque measurements (25j 
showed that a nematic transition occurs at a temperature 
much higher than the structural transition and extends 
to a region where the SDW phase is already far away 
but where superconductivity still occurs; consequently 
superconductivity is realized inside the nematic phase in 
a broad region of the phase diagram. 

In this paper, we study superconductivity due to or- 
bital nematic fluctuations, employing an effective two- 



band model which possesses an orbital nematic instabil- 
ity. We solve the linearized Eliashberg equations and 
take also the renormalization factor into account, which 
is often neglected in literature, but is essential to our 
study. We find that superconductivity is driven by 
long-wavelength orbital fluctuations and not by short- 
wavelength ones as discussed in Ref. [22|. 

Our model Hamiltonian has the form H = H + Hi 
where the interaction part H\ is given by 



Hi 



E 



Ui-Tli 



(1) 



The difference density operator is defined by nj_ = 
n-n - n l2 with the density operator n la = J^a c laa c ^aa- 
i and a are site and spin indices, respectively, and a — 
1,2 is a band index, g is a coupling constant which is 
considered as a parameter in our model. An expression 
for Hq suitable for pnictides is [26[ 



H = E E lf C L<r C k/3a, 
k,cr,a,/3 



(2) 



with e k * = —2tiCOsk x — 2t2COsk y — At^coskxcosky, 
e 2 . 2 = —2t2Cosk x — 2ti cos k y — At^coskxCosky, e k 2 = 
—4^4 sin k x si n ky . Reasonable values for the hopping am- 
plitudes are [if t = -h,t 2 /t = 1.5, t 3 /t = -1.2,t 4 /i = 
—0.95, which we will also use in our calculations; the 
chemical potential is fixed to be 0.6t. In the following 
energies are always given in units of t. 

The above model shows for negative values of g a ne- 
matic phase transition, i.e., a line which separates states 
where the expectation value is zero from those 

where it is nonzero. Figure 1 shows the result of a mean- 
field calculation of the phase diagram in the T-g plane. 
Disregarding superconductivity the region to the left of 
the solid line is occupied by the normal state whereas to 
the right side of this line a nematic state is realized. The 
nematic state is homogenous, i.e., is independent 



2 



of i, except for a small region bounded by the line with 
crosses and the dashed line in Fig. 1 where a modulated 
nematic state is found. 

In the following we will study superconductivity from 
nematic fluctuations. We consider the usual Fock di- 
agram for the electronic self-energy where the bosonic 
propagator describes nematic fluctuations. Since super- 
conductivity is a Fermi surface effect, we may restrict 
the momenta to the region near the Fermi line of each 
pocket. In our case there are two hole pockets near the 
F- and M-point and two electron pockets near the X- and 
Y-point. We denote them by i=l..A, respectively, see the 
inset in Fig. 1. Assuming that the superconducting or- 
der parameter has s-wave symmetry for each individual 
pocket the Eliashberg equations for the gap Ai(iu! n ) and 
the renormalization function Zi(iu} n ) read as 



j,n' 



iUJn') 



Zi(iu) n ) = 1 - 7r T Nj 



On' gijji 
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FIG. 1. (Color online) Phase diagram in the T-g plane for 
go = 0. The thick solid line separates the normal (N) from 
the orbital nematic (ON) phase which at low temperatures is 
modulated as indicated by the crosses. The dashed line is the 
instability line of the normal with respect to the homogeneous 
nematic state. The circles, smoothly joined by a thin solid 
line, separate the superconducting (SC) from the normal and 
nematic states. The inset shows the Fermi lines of the four 
pockets in the normal state. 



Lo n is a fermionic Matsubara frequency and Nj is the den- 
sity of states of pocket j at the Fermi energy. gij(iui n — 
iu) n i) is obtained from the microscopic pairing potential 
W a p(k., k', iv n ) by performing an average of k and k' over 
the Fermi lines of pockets i and j, respectively, and by 
summing over a and /3. W is given by 

W a/3 (k,k'>„) = ^(k,k')<?(k-k'> n ), (5) 

where v n denotes a bosonic Matsubara frequency. V rep- 
resents the form factor 



^(k,k') = £tft 7 (k)(r 3 ) 7S EMk'). 

7,5 



(6) 



U a p(\i) is a unitary matrix which diagonalizes Ho, t% a 
Pauli matrix and g(k, w n ) is given by 



g(k, iv n ) 



5 2 n(k, iv n ) 



9a- 



(7) 



- gFI(k,ii/, 

The first term in Eq. ([7]) is the retarded interaction me- 
diated by nematic fluctuations, the second one is an in- 
stantaneous term accounting for repulsive Coulomb inter- 
actions. II stands for a single bubble of non-interacting 
nematic particle-hole excitations. The transition tem- 
perature T c to superconductivity is obtained from the 
condition that the largest eigenvalue of the matrix 

M(i,n;j,n') = -irTNjgij(iuj n - iu> n ,)/\u> n , \/Zi(iu> n ) 

(8) 

is equal to one. 

The momentum-averaged pairing potentials gij(iv n ) 
are important ingredients in the calculation of T c . Fig- 
ure 2 shows the dependence of gij on the Matsubara fre- 
quency is n , treated as a continuous variable, for various 



temperatures. In the case of c/44 the momentum is av- 
eraged over the Fermi surface of pocket 4. Since the di- 
ameter of this pocket is small, only small momenta con- 
tribute to the average. The curves for g^iuJn) of Fig. 2 
(a) are always in the normal state. With decreasing tem- 
perature the correlation length of nematic fluctuations 
increases which implies that 344 (iuj n ) increases monoton- 
ically with decreasing temperature. At low temperatures 
g4i(iv n ) assumes huge values at v n = and then decays 
very fast with increasing frequency on an energy scale 
much smaller than t. In this frequency range 344(1^) 
varies only slowly with temperature which indicates that 
the attractive pairing interaction is rather insensitive to 
temperature and to critical fluctuations. For g = —1.8 
[Fig. 2 (b)], the nematic instability occurs at the temper- 
ature T n « 0.125. Here gii{w n ) first increases monoton- 
ically with decreasing temperature until T n is reached. 
Entering the nematic state, ga^Vn) is suppressed mono- 
tonically due to the development of nematic order. Ex- 
cept for this temperature dependence, g^{iv n ) is char- 
acterized by similar features as in Fig. 2 (a): large static 
values over a wide temperature region, a fast decay on 
an energy scale much smaller than t, and a weak temper- 
ature dependence at high frequencies. 

The insets in Fig. 2 show the dependence of the non- 
diagonal interaction 314 on v n . In this case the momen- 
tum average involves large momenta in the effective in- 
teraction and the spectral function of the nematic fluc- 
tuations extends over a large frequency region of several 
i's. As a result the decay of 514 with v n is much slower 
than in ^44 and characterized by an energy scale of t. 



(b) g=-l.8, g =0 




FIG. 2. (Color online) Retarded pairing interactions ga 
(main diagrams) and gi4 (insets) as a function of v n for dif- 
ferent temperatures for g = —1.7 (a) and g = —1.8 (b). In 
the case (b) the nematic state occurs below T — T n ~ 0.125. 

On the other hand the low frequency values of gu are 
much smaller than those of 344. In contrast to the case 
of 344, the temperature dependence of 1714 is very weak. 
We also found that the curves for 344 and gu in Fig. 2 
are representative for all diagonal i—j and non-diagonal 
i^j interactions, respectively. The only exception is 334 
which becomes practically zero due to the matrix element 
14/j(k,k') inEq. ©. 

We have solved Eqs. ©-([H numerically using up to 
2546 Matsubara frequencies to get convergence. Fig- 
ure 3 (a) shows the largest eigenvalue of the matrix 
M(i,n;j,n') as a function of temperature for go = 
and g = —1.7. Keeping only the diagonal 11,22,33,44 
(non-diagonal 13,14,23,24 or 12,34) contributions in this 
matrix, one obtains the dotted (short-dashed or long- 
dashed) curves; the solid line contains all contributions. 
The two curves corresponding to non-diagonal contribu- 
tions are small down to low temperatures. They never 
cross the value one in the chosen temperature interval 
and thus do not lead to superconductivity. The closeness 
of the solid and dotted lines proves that the total diag- 
onal contribution is much larger than that of the non- 
diagonal terms. T c is thus determined to a large extent 
by excitations with small momentum transfer. 

To understand the connection between T c and the fre- 
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FIG. 3. Temperature dependence of the largest eigenvalue 
using only diagonal (dotted line), non-diagonal (short and 
long dashed lines) or all pairing interactions (solid line). 

quency dependence of the pairing potentials g^j (Fig. 2), 
we first note that the zero frequency value of the pair- 
ing potentials with i—j does not enter the largest eigen- 
value when it crosses one. This follows immediately 
from Eqs. ©-Q, highlighting the crucial role of Zi(iui n ). 
Thus the strong attractive part in the pairing potentials 
between and 27rT c does not contribute at all to the 
instability towards superconductivity. At large temper- 
atures mainly the flat high-frequency tail of gij matters 
and already here the diagonal contribution accounts for 
most of the largest eigenvalue. Decreasing the temper- 
ature the low-frequency behavior of gij becomes more 
and more important because of the substantial increase 
of the number of Matsubara frequencies contributing to 
the right-hand side of Eq. <|3j> . This is the case also for 
Fig. 2 (b) although the low energy weight is suppressed 
at lower temperature in the nematic phase. As a result 
the huge increase of attraction towards low frequencies 
in the diagonal contributions dominates the total pairing 
potential, pushes the largest eigenvalue up to one and 
thus determines T c . 

Figure 3 (b) shows the largest eigenvalue for g = —1.8 
where the nematic phase is stable below T n ps 0.125. The 
non-diagonal contributions are again small and do not 
cross the value one. The diagonal and also the total con- 
tribution behave similar as in Fig. 3 (a) except in a tern- 
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perature region around T„ where the curves approach one 
from below. This findings follows directly from Eqs. ([3])- 
(U). At the transition point T n , ga(iv n — 0) and Zi(iuio) 
tend to minus and plus infinity for i — 3 and 4, respec- 
tively; wo is the lowest fermionic Matsubara frequency. 
This implies that the largest eigenvalue approaches one 
from below, but never crosses one. This could be in- 
terpreted in terms of enhanced superconducting fluctu- 
ations near T n . However, the assumptions underlying 
Eqs. ©-(JI} could break down near T n , and non-Fermi 
liquid features might appear similar as in Ref. I27I 

The pairing potentials gij{iv n ) are for go = attrac- 
tive for all pocket indices i,j and all Matsubara frequen- 
cies. As a consequence the components of the eigenvec- 
tor Ai(iv n ) belonging to the largest eigenvalue have the 
same sign for all pockets and frequencies. This means 
that for go — only the s ++ symmetry for the order pa- 
rameter can be realized. It is interesting that the second 
largest eigenvalue lies only a little below the solid line in 
Fig. 3 in the normal state. The corresponding eigenvec- 
tor has d-wave symmetry, i.e., Ai(iuj n ) = A2(iui n ) = 
and Asiiuin) = — A^iujn) for all frequencies. This result 
is plausible because g^ = 344 dominates the total pair- 
ing potential and #34 is very small. The pockets 3 and 4 
therefore decouple approximately which leads to a near 
degeneracy of the s- and the d-wave states. 

The filled circles in Fig. 1, smoothly joined by a thin 
solid line, show T c as a function of g. With increas- 
ing —g, T c monotonically increases because the nematic 
transition at g n w —1.72 is approached and the nematic 
fluctuations increase. T c decreases somewhat when pass- 
ing through the transition to the nematic state. Increas- 
ing —g further enhances T c monotonically. Although 
the fluctuations are suppressed deep inside the nematic 
phase, they are sufficiently strong to drive a supercon- 
ducting instability, leading to a coexistence of the super- 
conducting and nematic states at lower temperatures. In 
the nematic phase, the Fermi pocket 3 (4) can disappear 
for (nj_) > (< 0) at low temperatures, which is indeed 
the case for — g > 1.8. However, T c is not affected much 
by such a Fermi surface reconstruction. 

Zi(iui n ) assumes in general values much larger than 
one at low frequencies, which indicates that the strong- 
coupling limit of superconductivity applies. Taking the 
free value Zi(iuj n ) = 1 in Eq. © would lead in general 
to unphysically high critical temperatures and thus to 
incorrect conclusions. Approximating the frequency de- 
pendence of ^44 by a two-square- well model [11] the cutoff 
frequency would be about w c ~ 0.05. Since on the aver- 
age T c ~ 0.03 we obtain the ratio T c /lj c ~ 0.6 which is no 
longer small compared to one and thus again points to the 
presence of the strong-coupling limit. T c = 0.03 corre- 
sponds to about 50K for the hopping t = 150 meV. Since 
T c is determined by low-frequency, intra-pocket nematic 
fluctuations with an energy scale of uj c , which is much 
smaller than characteristic electronic energies, Migdal's 




FIG. 4. Dependence of T c on the repulsive, instantaneous 
interaction go- 



theorem should be rather well satisfied. This justifies our 
neglect of vertex corrections in Eqs. ©-(HI). 

In Fig. 4, we show the dependence of T c on the instan- 
taneous repulsive interaction go, which will be present in 
real materials and might hinder nematic superconductiv- 
ity. For the case g — —1.7 superconductivity evolves from 
the normal state and T c is only moderately suppressed 
by go- The components of the eigenvector belonging to 
the largest eigenvalue have at low frequency the same 
sign, hence s++ symmetry [29J, but change the sign at 
a finite frequency. This sign change leads to a partial 
cancellation of the repulsive terms in Eq. ([3]), which al- 
lows a non- vanishing solution for Ai(iu> n ) even for a large 
go. For the case of g = — 1.9 in Fig. 4, superconductiv- 
ity arises inside the nematic phase. Comparing the two 
panels in Fig. 2 indicates that retardation effects are in 
the nematic phase weaker than in the normal phase at 
low temperatures. The cancellation of strong repulsive 
terms in Eq. ([3|) by sign changes of the eigenvector is 
then impeded which explains the faster decay of T c and 
the vanishing of T c with increasing go- 

In iron pnictides, nematic fluctuations seem to occur 
even far away from the SDW phase [25|. The present 
theory is applicable directly to such a case. While s ++ 
symmetry is always the leading instability, it is nearly de- 
generate with a d-wave state for the instability from the 
normal state. Hence additional small interactions such as 
residual spin fluctuations may be sufficient to change the 
energetics. Upon approaching the SDW phase, it is likely 
that spin fluctuations also play a role for superconductiv- 
ity 3CJ, [31 1 . On the other hand, recent Raman scattering 
measurements revealed directly nematic fluctuations 
close to the SDW phase. One therefore expects quite 
generally a competition between spin and nematic fluc- 
tuations. The nematicity may thus play an important 
role to understand high-temperature superconductivity 
observed in iron pnictides. 
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